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Abstract 

We investigate the limiting distribution of the fluctuations of the max- 
imal summand in a random partition of a large integer with respect to a 
multiplicative statistics. We show that for a big family of Gibbs measures 
on partitions (so called generalized Bose-Einstein statistics) this distri- 
bution is the well-known Gumbel distribution which usually appears in 
the context of indepedent random variables. In particular, it means that 
the (properly rescaled) maximal energy of an individual particle in the 
grand canonical ensemble of the d-dimensional quantum ideal gas has the 
Gumbel distribution in the limit. 

We also apply our result to find the fluctuations of the height of a 
random 3D Young diagram (plane partition) and investigate the order 
statistics of random partitions under generalized Bose-Einstein statistics. 

1 Introduction 

The link between the energy distribution in grand canonical ensembles of a 
quantum ideal gas and the asymptotic behavior of random partitions of integers 
was established in several papers. For example, the problem about the limit 
shape of the energy distribution in such ensemble was considered in ^1 . In 
this paper we are interested in a more subtle question — the limit distribution of 
the maximal energy of a particle in the grand canonical ensemble of the ideal 
Bose gas of dimension d with the Gibbs statistics. We give the answer (see 
Theorem 13 below) , but the question is involved in a general context which has 
some history. 

Consider the set of all partitions V = U7'(n) of all natural numbers n with 
the so called multiplicative (Gibbs) statistics, — it means that, firstly, occupation 
numbers are independent with respect to this statistics and, secondly, restric- 
tions of the statistics on the subset of all partitions of the given number n does 
not depend on the parameter of the Gibbs measure (activity). We give precise 
definition below. The simplest example of a multiplicative statistics is the so 



called Poissonization of the uniform measures on partitions of integers; a more 
general example is the Gibbs measures of the quantum ideal gas, either Fermi 
or Bose (see (EI)- It is possible (under some special conditions) to find the 
limit shape of partition, i.e. the distribution of energy |14j . The next problem 
is to consider the distribution of the fluctuations. In the bulk of summands the 
fluctuations have a Gaussian distribution (see |15[I18) '). But the fluctuations on 
the "ends" of partition, i.e. the fluctuations of the maximal summands, have 
completely different form. We can compare this situation to the same question 
for Plancherel measure where the distribution of the fluctuations of the longest 
increasing subsequences (or the maximal row of a random Young diagram) with 
respect to the Plancherel measure is the Tracy- Widom distribution (0). For 
our problem we get a so called Gumhel distribution with the distribution func- 
tion e~'^ which occurs also in the theory of extreme values of independent 
random variables (see, e.g., 0), but its appearance in our case has a different 
nature. 

First time this distribution appeared in similar problems in the pioneering 
work of Erdos and Lehner ||6j where they proved that the maximal summand in 
a random uniformly distributed partition of a large positive integer n is approx- 
imately (\/6nlogn)/(27r) and found a right scaling and a limiting distribution 
for the maximal summand m(A) in a random partition A if the partition is taken 
uniformly among all partitions of n. Exactly, their result is that for all a; g R 

1 r m(A) - (v/^logn)/(2c) -i -^-'e-"- 

where v(n) is the number of all partitions of n. The proof presented in |H] is 
based on some combinatorial estimates for numbers of partitions with certain 
properties, asymptotic relations based on these estimates and the inclusion- 
exclusion principle. 

Similar behavior is known for other measures on partitions, for example, 
for the uniform measure on set partitions. Namely, take at random a partition 
of the set {1, . . . ,n} and look only on block sizes, which give a (nonuniform) 
distribution on partitions of an integer n. In this case the limiting distribution of 
the maximal block size is a discrete approximation to the Gumbel distribution, 
see ini. Section IV. 5. 

In this note we present a different approach to this problem allowing us to 
generalize this result to a more general family of measures on partitions. We use 
Poissonized measures and prove all our results only for them. These results for 
Poissonized measures can be extended to the corresponding statements for orig- 
inal measures. General principle allowing to do it will be presented in another 
paper. 

Now we give the precise description of what is called multiplicative statistics. 
Let Tin) be the set of partitions of an integer n and V — U„>o7^(?^) be the set 
of all partitions, as above. For A — (Ai, . . . , A^) G V, \\ > ■ ■ ■ > Xe > 0, we 
write rfe(A) = {i : Xi = k} for the number of summands k in partition A, so 
called k-th occupation number. Consider a sequence of functions fk{z), > 1, 
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analytic in the open disk D = {z G C : |z| < R}, i? = 1 or i? = oo, such that 
/fc(0) — 1 and assume that (i) the Taylor series 



have all coefficients Sk{j) > and (ii) the infinite product 

oo 

F{x)^l[h{x') (2) 

fe=i 

converges in D. Then one can define a family of probability measures fJ-x, x G 
(0, R), on the set of all integer partitions V in the following way: put 

M.({Ae7':r,(A)=j})^f^^, (3) 

and assume that different occupation numbers are independent. Note that in or- 
der to specify measure iix it suffices to fix F{x) along with its decomposition 
At the same time just specifying F{x) is not enough. 

Definition. The family of measures fXx defined by (|3Jl and satisfying condi- 
tions (j) and (ii) is called a multiplicative family of measures (or a multiplicative 
statistics) on partitions. In this case we say that the statistics is determined by 
decomposition 

This notion in a general context was introduced in J^l; similar technique 
was exploited earlier by Fristedt [7] for the uniform measure on partitions. The 
key feature of these measures is the fact that conditional probability measures 
induced on ■p(n) do not depend on x, for all n. We denote these conditional mea- 
sures by /x" = Hxl^f^^y Thus measures fix can be considered as a Poissonization 
of measures /x" which makes random variables independent. Many state- 
ments which hold for fix hold for fj," too, see, for instance, |13l Hj: at the same 
time, the investigation of measures fix is simpler because are independent. 

We call the family {V{n), fj,") small canonical ensemble of partitions^ and the 
family (VyfJ-x) — grand canonical ensemble of partitions, in view of similarities 
with statistical physics, see [Ol. 

We restrict our attention to the special case of a so-called generalized Bose- 
Einstein statistics defined by Q with fk{z) = 1/(1 ^ z)^'' . Thus equation ((2Jl 
specializes to 

^^(-) - n 713^' (4) 

fc=l ^ ' 

b — More exactly, first we treat only the case when 6^. = ck^ for some 

c > and /3 > — 1. In this case direct calculations can be made to show that 

^The notion microcanonical ensemble is usually used for the set of partitions n = Ai -|- ■ • ■ + 
Xl with both weight n and length £ fixed, so we introduce the term small canonical ensemble 
for an ensemble with fixed energy but varying number of particles. 
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the limiting distribution of a maximal summand is the Gumbel distribution, 
as done in Theorem ^ Note that the Poissonization of the uniform measure 
on partitions belongs to this family (all bk = 1 in this case). For the uniform 
measure the distribution of a maximal summand coincides with the distribution 
of the length of a partition (that is the number of summands) , but in a general 
case these distributions are essentially different. 

For integer bk these partitions are often called colored partitions because the 
measure on the small canonical ensemble is induced by the uniform measure on 
partitions with additional structure: each summand k can be colored in one of bk 
colors, and two partitions are identical when number of summands of each size 
and color coincide, see pp. However while the requirement that bk are integers 
is natural from the combinatorial point of view, it is not needed analytically 
and we consider all real positive bk- 

Our main result is presented in Section|21 It concerns the distribution of the 
maximal energy of an individual particle in the large canonical ensemble of the 
d-dimensional quantum ideal gas (see, e.g., for background). We show 
that after a suitable scaling the limiting distribution of the maximal energy of an 
individual particle is the Gumbel distribution. The exact statement is given in 
Theorem|21 Note that while the grand canonical ensemble of the quantum ideal 
gas is the partial case of the generalized Bose-Einstein statistics. Theorem ^ 
can not be directly applied since bk do not have a form ck^ even asymptotically. 

In Section ^ we sketch the application of our result to another example, 
namely to the the distribution of the height of 3D Young diagrams. As above, 
we consider only the Gibbs measure on these objects. In Section|3we conjecture 
the limiting behavior in small canonical ensembles of partitions. The last section 
is devoted to the investigation of the order statistics, or, in other words, of the 
spacing between first largest summands in partition. 



2 Maximal summand 

We consider the special case of multiplicative statistics, namely measures /ib^^, 
defined by ® with/fc(z) = l/{l-z)'"', b = {bk} = {ckl^}k>i with /3 > -1. The 
convergence radius of fk and F is R = I. We denote the maximal summand in 
a partition A = (Ai, . . . , A^), Ai > • • • > > 0, by to(A) = Ai. 

Theorem 1. Let the measure ^b^x be defined by Q where 

k>i k>i ^ ' 

for some c > and /? > — 1 . Then, for all t (^M 

lim Hb.xiX e V ■■ to(A)(1 -x)- A{x) < t] = e^^'\ (5) 

where 

A{x) ^ (/3+l)|log(l-x)| +/31og|log(l-x)| +/31og(/3+l) + logc. (6) 
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Remark. The standard form of a limit theorem is to find Um/i^rjA G V : 
"^^^b{x)^^^ — Expression looks differently but can be rewritten in this 
form taking a{x) — A{x)/{1 — x) and b{x) — 1/(1 — x). 

Proof. We shall not just verify the statement of the Theorem but also show how 
it can be deduced. First, note that the probability fib,x{^ G "P : m(A) < M} 
tends to for fixed M as x — > i?, so in order to get sensible results we should 
take M depending on x. More exactly, in order to get the limit theorem in 
form ||SJ| we take 

M = M{x, t) = a{x) + b{x)t, (7) 

where o and h are functions of x presumably growing to infinity as a; ^ 1 , and t 
is a parameter. Since we want M{x, t) grow to infinity for all fixed t, it follows 
that h{x) = o{a{x)) as a; ^ 1. 

Since measures ^b,x are more adjusted to work with occupation numbers r^, 
we reformulate the question in terms of r^. in the following way: for any M > 1 

Mb,.{A G r : mix) < M} = //,{A : r^iX) = for fc > M} = — (8) 

k>M ■'''^^ > 

Thus we just have to find functions a(-) and 6(-) such that, as x ^ i?, the 
product of l/fk{x'') taken for k > M{x,t) = a{x) + b{x)t tends to some non- 
degenerate distribution function. 

We take minus logarithm of (jSJ to conclude that 

-\oginbA^:m{X)<M{xA)}) = - ^ 6fc log(l - x'^). 

k>M(x,t) 

Let us assume that 

Af(a;,t) >p|log(l-a;)|/|logx| (9) 

for some p > 0; under this assumption we have < supj,>„(^_j-) x'' < C{1 — x)p 
for some C > 0. So, using ^ and putting bk = ck^ we can calculate the sum 
above explicitly: 



2k 



- log(^5,^{A : m(A) < M{x,t)}) = ^ ck'^x'' + o( ^ kl^x 

k>M{x,t) k>M{x,t) 

= «^(1 + 0(1/M(M))) (10) 

as x -^^ 1 with fixed t. We take a logarithm once more and arrive to 

log(- log(/ib,,{A : m(A) < M{x,t)})) 
^-M{x,t)\\ogx\ + \ log(l - x) I + /3 log M {x,t) + \ogc + 0{l/M{x,t)) (11) 

as x ^ 1. 
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We are seeking M{x,t) in the form {Tj), and taking b{x) = l/|logx| in Q 
we can rewrite as 

log(-log(Aib,,{A : m(A) < M{x,t)})) 

= -a(a;)|loga;| + | log(l - a;)| + pioga{x) + logo - t + 0{l/a{x)) . 

So if we could find a(x) such that 

-a{x) I log xl + I log(l -x)\+ (3 \oga{x) + log c = (12) 

then fj,b.x{^ ■ m{X) < AI{x,t)} would tend to * for fixed i as a; — > 1. We 
are searching the solution of (|12|l in a form 

(/3+l)|log(l-x)|+ai(x) 

a a; = n \ 

|loga;| 

where ai{x) = o(|log(l — a;)|). After a substitution into (|12|) we immediately 
see that ai{x) = /Slog | log(l — a;)| + logc + o(l). It remains to use the relation 
I loga;| '--^l — xasx— >lto get (jSJ. 

Now note that for this choice of a(-) and &(•), the assumption is satisfied 
iff /3 > —1, which justifies all computations above. □ 



3 The quantum ideal gas 

Here we apply our results to obtain the limiting distribution of the maximal 
energy of an individual particle in the quantum ideal gas in M'*. (The most in- 
teresting case is of course d — 3 but the method works for all d > 1). We describe 
briefly the connection of the quantum ideal gas to the partition theory; for a de- 
tailed exposition see, e.g., [HI El- To each configuration in the phase space there 
corresponds a partition of an integer n: in the suitable units, the energies of indi- 
vidual particles become summands in the partition and the energy of the whole 
system becomes n. There are jd{k) = #{(fci, . . . , kd) & : kf + ■ ■ ■ + k"^ = k} 
distinct positions of particles in the phase space such that a particle in these 
positions has (rescaled) energy k. Consequently, the measure on partitions in- 
duced by the grand canonical Gibbs measure on the quantum d-dimensional 
ideal bosonic gas is a multiplicative measure determined by the decomposi- 
tion (0J with bk — jd{k). This observation was used by Vershik ^|E] to find a 
distribution of energy among particles in these settings. The following statement 
describes the behavior of the maximal energy of an individual particle. 

Theorem 2. Let measure fi^.x be the multiplicative measure on partitions orig- 
inating from the quantum ideal Bose gas, i.e. it is defined by (PJ where 

F(.)^n^(^^)-n n_L.w - 

k>l k>l ^ ' 
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Then, for all t G 

lim /i^{A e V : m(A)(l - x) - A{x) < t} ^ e'"'' 

where 

A{x) = f |log(l -x)\+^ log|log(l - x)! + I log I + log r(j','^,) . 

Classical results on jdik) state that jd(fc) = 0{k'^^^^^) for d 7^ 4 and j4(fc) = 
O(fclogfc), see jS]. Nevertheless, jd{k)k^''^^'^ has no limit. Thus, Thcorcmncan 
not be directly applied to obtain results on limiting behavior of the maximal 
energy of an individual particle. But the ideas used in its proof still work. 

Proof of Theorem O We are going to get estimate analogous to (fTIHl used in the 
proof of Theorem n which was crucial to get that result. Once we have such 
estimate the rest of the proof almost literally repeats the proof of Theorem ^ 

Take M{x,t) = \[A{x) + loga;|] ; then supj.>j,^(^ x'^ 0. As in the 
proof of Theorem n we obtain 



-log(/Xb,,{A:m(A)<M(a;,i)})= ^ ]d{k)x^+o(^ ^ jd(fc)a 

k>M{x,t) k>M(x,t) 

(13) 

Both sums above can be treated in a similar way. We show that the first sum 
tends to as x — > 1; replacing a; by in suitable places shows that 0{-) 

vanishes. 

Denote Jd{k) — J2i<kjd{i)- It is well known (see |H| E|) that Jd{k) = 
Cdk'^^'^ + Ed{k) where Cdk'^^'^ — r(d/2+i) ^''^^ volume of a d-dimensional 

ball of radius Vfc and the error term Ed{k) = 0{k°''') where ad can be taken as 
follows: ai = 0, a2 = 1/3, as = 3/4, a4 = I + S {6 > 0) and ad = d/2 — 1 for 
d > 5 (better estimates are known but we do not need them). 

For any M > 1 we can write 

)x'' 



]d{k)x'^ = ('^'i(^) - '^dik ' 1)) 

k>M k>M 

= CdY {k"'^ - - 1)'^')^' + {l-x)Y Edik)x'' ~ Ed{M - 1 



X 

k>M k>M 

^dx^^M(d-^y^ ^ 

= Crf- + Rd{M), 

2 |loga;| 

where the error term Rd{M) can be estimated as follows: \Rd{M)\ < |i?^(Af)| - 
\RdiM)\ with 

R'diM) ^CdY {k^'^ -{k- iY'^)x^ Cd\ ;! , 

k>M ' ' 



R'lliM) = {l^x)Y Ed{k)x'' - EdiM - 1 



k>M 



7 



Using the integral approximation for the sum one can find that R'^{M) < 
Kix'^^M^'^-'^''/^ and estimates on Ed show that i?Jj'(M) < K2X^ M"'" for some 
Ki, K2 > 0. Taking M = M{x, t) = [{A{x) + t)l\ logx|] we see that the lead- 
ing error term is R'^{M) and that for fixed t the inequality \Rd{M{x,t))\ < 
K{1 - xY''^-'^'' holds for some K > Q. It remains to check that our choice of 
M{x,t) imphes that 

^d x^^^'^^^){M{x^t)f-''" 

2 |10gx| 

Almost the same calculation verifies that the argument of O(-) in 1)13(1 vanishes 
as a; ^ 1. □ 

One can also consider the Fermi-Dirac rf-dimcnsional quantum ideal gas, 
which also induces the multiplicative statistics on partitions with fk{z) — (1 + 
zy'i'^^\ In the grand canonical ensemble it has exactly the same limiting be- 
havior of the maximal energy of a particle. It can be seen analytically since 
the limiting behavior of the maximal energy depends only on the first Taylor 
coefficients of functions fk, and they are the same for the Fermi-Dirac and 
Bose-Einstein statistics. It worse noticing that most of statistical properties of 
Fermi and Bose gases are different. 



4 A section of 3D Young diagram 

Let us consider a set of 3D Young diagrams (or plane partitions) of weight TV. A 
3D Young diagram is a Z-|_-valued function h{u,v) of two arguments u, u e 
such that it has a finite support, is non-increasing in both arguments and if it 
is discontinuous in a point (u, v) then either m or w is integer. A weight of a 
diagram h{u, v) is 

N{h{-,-)) = // h{u,v)dudv 



and it is obviously an integer. The graph of h{u,v) in M'^ is an upper bound 
of a set which can be constructed from unit cubes in the same way as ordinary 
2D Young diagram is constructed from unit squares (or boxes); this analogy 
explains the name of these objects. We denote the set of all 3D Young diagrams 
by V^iD- 

Given a number x e (0, 1) one can consider a probability measure on the set 
of all 3D Young diagrams with probability of any diagram /i(-, ■) is proportional 
to x^^'^'--'''^^ The generating function for numbers psiN) of 3D Young diagrams 
of weight N is well known: 

E ^3(^)-^ = n 



Af>0 k> 



^r(l 



see, e.g., The existence of a limit shape for 3D Young diagrams was 

proved by the first author; the exact formulas were found later by Kenyon 
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and Cerf [3] using a variational technique and further investigated by Okounkov 
and Reshetikhin JI] . Thus any section of a graph of 3D Young diagram has a 
hmit shape also. In particular, a section by the plane u — v can be considered 
as the Young diagram of an ordinary partition (up to a factor ^/2). It turns 
out that a distribution of these diagram will be exactly the distribution defined 
by l@J with hk = k (we should treat h axis in the 3D case as t axis in 2D). The 
correspondence between colored partitions and diagonal sections of 3D Young 
diagrams is rather complicated and includes Bender-Knuth bijection between 
random infinite integer matrices and pairs of semi-standard Young tableaux and 
the correspondence between these pairs and 3D diagrams. It was introduced by 
the first author in his talk jJSj, its detailed description and consequences will 
be presented in a separate paper. 

Thus, direct application of Theorem pleads to the following result: 

Theorem 3. Let fi^ be the Gibbs probability measure on Vsd, the measure 
of a particular diagram h{-, •) is proportional to x^^'^^'''^\ Then, for all t G M, 

limfiJhi-,-) eVsD ■■ h{0,0){l-x)-Aix) <t} -e"^"*, (14) 

where 

A{x) = 2|log(l - x) \ +log|log(l - x) \ -flog 2. (15) 

5 Statements for small canonical ensemble 

The corresponding results in small canonical ensemble of partitions, that is 
results for measures can be obtained by taking x = x{n) depending on n so 
that the expected weight of partition {J2k ^^fe) ^'-'^ measures considered 
in Theorem ^ it can be achieved by taking 

x = xin)^l-^^{[^^^±M^ 
V n 

(C is the Riemann zcta function), see ^S]- We say that grand canonical and 
small canonical ensembles are equivalent for some functional G on partitions if 
the distributions of the functional G w.r.t. /i" and w.r.t. ^ix{n) are asymptotically 
the same. It seems that the ensembles are equivalent for the functional of 
rescaled maximal summand in partition. We shall return to this question in 
another paper. 

In the assumption that ensembles are equivalent Theorem^yields the follow- 
ing result: if measures fi"' are conditional probability measures induced on V{n) 
by the generalized Bose-Einstein measures fix determined by decomposition (0J 
with bk = ck^ , then 

lim ^"|a e Vin) : ^W ^r(/3 + 2)C(/3 + 2) _ < ,1 ^^g) 

n-*oo L V n J 
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where 



+ , 

■ log n + p log log n 



/3 + 2 



The terms in the second line above do not depend on n so they constitute a 
constant correction term, while the terms in the first line show how maximal 
summand is growing with the growth of n. 

In two particular examples considered above, the exact computation can be 
made. For the case of the quantum ideal gas of total energy n (in the suitable 
units so that n is integer), we should take 



x = x{n)^l- 

in the grand canonical ensemble to get the best approximation of the small 
canonical measure /i". Thus, under the assumption of equivalence of ensembles, 
in the small canonical ensemble of d- dimensional quantum ideal Bose gas 

^lim m"{a G V{n) : m{X)[^^^^^^^^y^^'^'^ - A{n) < t} = e^"* 

with 

^" = ■5|2logn+ ^loglogn 

+ 2(lW i°g I + ^ i°g - ^2 iog(C(rf/2 + l)M - logr(d/2 + 1). 

Similarly, it follows from H16() that a height of a typical 3D Young diagram 
(i.e. h{0, 0)) of weight N behaves as 3(^^3)1/3 N^^'^ logiV as N ~* 00. More exactly, 
under the assumption of equivalence of ensembles, denoting by the uniform 
measure on 3D Young diagrams of weight N , we have 

^lirn^ v''[h: Qi^yfflh{Q, 0)-f log TV-log log iV-log |+| log2C(3) < - e"^" 



6 Order statistics 

Now let us consider the upper order statistics of random partitions, that is the 
sequence of the first d largest summands toi(A) > • • • > md{X) in a random 
partition A. We consider only the settings of Theorem ^ generalizations to 
other measures considered above are straightforward. 

Theorem 4. Let measures ^ be defined by with fk{z) = 1/(1 — z)^'' , b = 
{bk} = {ck^}k>i and (3 > —1. Let A(x) be defined by (Q. Then the distribution 
of rescaled upper order statistics (mi(A)(l — x) — A{x))_.^-^ ^ converges weakly 

to the distribution on M'* with the joint density exp(— e^*"* — X^iLi ''0 f^^ ^1 
■ ■ ■ > td and zero otherwise. 
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Proof. Let us fix > f2 > • • • > and let (5 > be such that intervals 
Ai = [ti, ti + S] are disjoint. Then the probability that each mi(A) after rescaling 
gets in the interval is 

/ix{A e V : m,(A)(l - x) - A{x) G A„ i = 1, . . . , d} = „ ^'^^ , (17) 
where 

(Recall that Sfc(l) = fcfe = cfc'^ and Sk{l)x'' / fk{x'') is the probability that rt = 
1.) Note that we estimated the denominator in the RHS of 1)17(1 while we proved 
Theoremn(see equations ((TT|) and and that it tends to 1/ exp(— e"*"*) as 

a: ^ 1. Let us estimate the numerator, that is sums Si. In view of HlOfl. 

= ^{{Mix,u fx'^''-^'''^ - {Mix,U + 5)fx'''^^'''+'^) 

where M{x, t) = + j^t. This choice of M{x, t) imphes that Si e^*' (1 - 
e"'') as X ^ 1 for fixed ti. Taking 5^0 limit finishes the proof. □ 

Remark. The same behavior of upper order statistics is known for samples of n 
i.i.d. random variables lying in the attraction domain of Gumbel distribution, 
see, e.g., 0. 
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